Conductivity between Luttinger liquids: coupled chains and bilayer graphene 
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The conductivity properties between Luttinger liquids are analyzed by exact Renormalization 
Group methods. We prove that in a two chain system or in a model of bilayer graphene, described 
by two coupled fermionic honeycomb lattices interacting with a gauge field, the transverse optical 
conductivity at finite temperature is anomalous and decreasing together with the frequency as a 
power law with Luttinger liquid exponent. 
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INTRODUCTION 



The many body interaction in fermionic systems can 
destroy the electron-like nature of the elementary ex- 
citations, a fact which can have deep consequences on 
the transport properties. This was pointed out first by 
Anderson wn0 g°t evidence that the coherent trans- 
port between Luttinger liquids due to particle hopping 
is strongly depressed with respect to the non interacting 
case. One of the first applications of this idea was, see 
@j 1 a- n explanation of the c-axis anomalous conductiv- 
ity between planes in high T c superconductors, with the 
assumption of Luttinger liquid behavior in the planes de- 
scribed by 2D square lattice Hubbard models; while in- 
triguing, this theory suffers the fact that no convincing 
theoretical evidence has been found up to now to sub- 
stantiate such assumption. Subsequently, the attention 
was focused to one dimensional systems, where Luttinger 
liquid behavior is surely present; in addition to a theo- 
retical interest, the main physical motivation was that 
coupled fermionic chains well describe quasi-one dimen- 
sional organic conductors, see e.g. At zero tempera- 
ture Renormalization Group or bosonization analysis ap- 
parently indicate that the hopping can destroyLuttinger 
liquid behavior in several regimes, see e.g. [5l4l7|. On 
the other hand at higher temperatures the system still 
shows Luttinger liquid properties l3-21|- The transverse 
conductivity between weakly coupled fermionic chains at 
temperatures and frequencies greater than the hopping 
was shown to be a (uS) ~ Lu a with a related to the Lut- 
tinger liquid exponent r); in particular the value of the 
exponent a was claimed to be a — 2r], in 1221 . where a 
tunnelling approach was followed, and in 12311 . by using 
the Kubo formula; on the other hand, in [24j the dif- 
ferent value a = 2r/ — 1 was found, by using dynamical 
mean field theory. The reason of these discrepancies re- 
lies on the fact that the computations were done in the 
so called Tomonaga-Luttinger or g-ology approximation, 
in which the fermions close to the two Fermi points are 
described in terms of massless Dirac particles. This sim- 
plifies the computations and allows the use of powerful 
techniques like bosonization, but introduces spurious ul- 
traviolet divergences in the conductivity which therefore 



needs a regularization, and different regularizations pro- 
duces different results (see e.g. the discussion in [24J after 
eq.(3)). 

In conclusion, there are at present no firm results on 
the conductivity between higher dimensional Luttinger 
liquids, and even in one dimension there are still am- 
biguities due to approximations and regularizations. In 
this paper the conductivity properties between Luttinger 
liquids at finite temperature are analyzed by the exact 
Renormalization Group methods developed starting from 
pBj : indeed the suggestion of using such techniques for 
this problem dates back to Anderson himself |l| but their 
full development required a long time. In the case of cou- 
pled spinless chains we get an exact expression for the 
transverse conductivity (the lattice furnishes the natu- 
ral cut-off), which is, at temperatures and frequencies 
greater than the hopping 



"(w n ) 



(i) 



if 



= 2 f(n 



i). The analysis is based on the im- 
plementation of Ward Identities in the Renormalization 
Group, with rigorous bounds for the corrections due to 
the lattice 26]. Note also that the computation of the 
parallel conductivity gives in this regime a^(ui n ) ~ w^ 1 , 
that is no anomalous exponent appears in the frequency 
dependence in that case. 

In addition to justifying the anomalous exponent pre- 
dicted by mean field theory in [24[ for the two chain prob- 
lem, the exact Renormalization Group methods can pro- 
vide for the first time evidence for anomalous transverse 
conductivity between bidimensional Luttinger liquids, for 
which bosonization cannot be applied. Electrons on the 
honeycomb lattice interacting with an U(l) gauge field, 
representing the retarded e.m. interaction or the effect 
of disorder or ripples [27j ]. have Luttinger liquid behav- 
ior. Indeed this system was first analyzed in the contin- 
uum Dirac approximation in the early work [28j], where 
evidence of Luttinger liquid was found based on second 



order perturbation theory. Later on, in 28, 30] Luttinger 
liquid behavior was established at any order and taking 
rigorously into account the honeycomb lattice, by imple- 
menting lattice Ward Identities in the Renormalization 
Group scheme. The behavior of the two point function 
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is similar to the one of the spinless chain; the wave func- 
tion renormalization has a power law with exponent rj. 
By coupling two interacting fermionic honeycomb lattices 
by an hopping term we get a model for bilayer graphene 
3l| . The zero temperature properties of such system are 
rather complex and still not completely understood, see 
e.g. 32H34j and the review 35(. However, as in the case 



and 



of coupled chains, at finite temperature and frequencies 
the Luttinger liquid behavior of the uncoupled system 
can reveal itself by the transverse conductivity; in partic- 
ular we will show that for temperatures and frequencies 
greater than the hopping 



(2) 



while is essentially constant in the non interacting case. 
This confirms for the first time in two dimensions the 
Anderson idea: coherent transport between Luttinger liq- 
uids is depressed with respect to non interacting systems 
or Fermi liquids. Moreover, the presence of an anomalous 
Luttinger liquid exponents in the frequency dependence 
of the transverse conductivity of bilayer graphene could 
be revealed in future experiments. 

The paper is organized in the following way. In §11 
we derive the transverse conductivity for the two chain 
model; in §111 we derive the transverse conductivity for 
bilayer graphene. In App A the computation of the con- 
ductivity in the non interacting case is presented. 

THE TWO CHAIN MODEL 

We consider a two chain model described as two one di- 
mensional interacting spinless fermionic systems coupled 
by an hopping term; the Hamiltonian is 



H = H 1 +H 2 + P 
where, for i = 1, 2, Hi = Hf> + V t 

L-l 
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(3) 



(0) 1 4 

Hi — ~ 9 ( a x + l,i a x,i + a x,i a x+l,i) 



x = l 
L-l 



Vi = -a ^2 v ( x ~ y) a t,i a x,i a t,i a v,i ( 4 ) 

x,y=l 



and 



L-l 



P = ~t X! ^ a t,l a xa + a x,2 a x : l\ 



(5) 



where a xi are fermionic operators and |«(a;)| < e - "' 21 '. 
Either A and t are assumed to be small. 

As usual we can introduce the interaction with an e.m. 
field with a Peierls substitution H — > H(A) with Vi{A) = 
V and, if A= (M,A X ), 

H?\A) = -J2 \{< + i/ A Ki+<^ iAla * + ^ (6) 



P(A) 



e' A * a. 



+ -iA A 
x.2 T" a x,2 e 



a x,l) ( 7 ) 



The parallel current is defined as 



dH(A) 
dAl 



j^+Aij^+0((A"y) 



(8) 



where j x 'f and j x f are called respectively paramagnetic 
and diamagnetic part of the current and are given by 



2 1 



i=l 

2 1 

2 ( a i+t,i a x,i a xd a x+l,i) 



The transverse current is defined as 



dH{A) 
dA± 



ti' X +A^ + 0((A^r) 



(9) 



(10) 



where j x '~ L and j x J,± are called respectively paramagnetic 
and diamagnetic part of the current and are given by 

t 



J 
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p.± 



D.± 



'X a x,l a x,2 a x,2 a x,l) 
t{ , a iA a x,2 +a t,2 a x.l) 



(11) 



Finally the fermionic density is p x = a^ x a\ tX + x d2,x- 
If p = (oj n ,p), ui n — ^S-n,p = -j;m, the transverse 
conductivity at finite temperature is given by 

0-/3 (Wn) 



lim [ 

u> n p->o 



where, if A = O x 



(tt x ;j^) + (j?> x )] (12) 

°~ gH Tpl)] 



Tr[e" 



with T 



'xi ' ' ' Xn , (A) = ~Tr[ e - f>H] 

being the time order product and T denotes truncation. 
An analogous definition holds for the parallel conductiv- 
ity. 

In order to compute the conductivity (112p it is conve- 
nient to introduce a Grassmann integral representation 
for the correlation; we introduce the following generating 
functional 

e W, A (^) = f p(#)e _V«0-B(A^) (13) 

where ^ x i are Grassmann variables, x = (xq,x), P(dtjj) 
is the Grassmann integration with propagator Si.jg(k) 

1 



5(k) = 



-iko + cos k — cospf 



(14) 



with ko = ^-{n + V is the interaction given by, if 
/ dx = /-a/2 dx ° 



2 r 



i=l 



dxip+M^t - t / dx(V'ilV' x ,2 + V'x,2V , x,l) 
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with u(x— y) = 5(x — y )v(x — y) and v is a counterterm 
which is introduced to take into account the renormaliza- 
tion of the chemical potential; moreover 



Defining Ht (x) the Fourier transform of ( ; ) we can 



write 



H t (*) = r 



Wt,A(A- L )| 



(16) 



The analysis of the functional integral (|T3|) will be done 
by Renormalization Group (RG), integrating smaller and 
smaller momentum scales. The hopping t introduces an 
intrinsic scale in the RG analysis. For scales greater than 
the (renormalized) hopping, it is natural to treat the hop- 
ping as a perturbation using the chain variables ^ , on 
the other hand at smaller scales the hopping cannot be 
considered a perturbation and it is convenient to use the 
variables 



&k,2 = $M 



V>k, 2 ] (17) 



in terms of which the free action is diagonal but the 
& k 1; 6 k 2 have different Fermi momentum. Note that the 
temperature acts an infrared cut-off so that for temper- 
atures not too small only the first regime is present. 

The first step of the RG analysis is the decomposition 
of the propagator g(k) as a sum of propagators supported 
close to the two Fermi points ±Pf and more and more 
singular in the infrared region, labelled by a quasi parti- 
cle index a = ± (labelling the Fermi points) and by an 
integer h < 0: 



g(k)=gW(k)+ J2 E^-Pf) 

h=h@ ex.— ± 

with p£ = (0,p F ), gi h) supported on < |k- p<* \ < 
2 h+1 and <?W(k) has support far from the Fermi points. 
Note that 2 h f ~ 7r//3; the fact that the temperature is 
finite implies that there is a finite number of scales. 

The RG integration procedure is defined recursively in 
the following way. Assume that we have integrated the 
fields ipf',ij}fl, we S et 



e Wt,\{A x ) = e F h (A x ) 



p( d ^(<h)^ e ~V (h \^^ h >)-B( h HA ± : V^*4>'~ h >) 



where P(d^- h ^) is the Grassmanian quadratic integra- 
tion with propagator given by 



ilex 



-iko + cospf — cos(fc + cxpf) 
(20) 



with Xh(M) is a smooth cut-off function with support 
|k — api?| < 2 h+1 and Zh is the wave function renormal- 
ization. The single scale propagator g£\x.) is obtained 
from g"a (k) replacing Xh(M) with /ft(k) with support 
2 h " 1 < |k'| < 2 h+1 , k = (fc , k'), k = k' + ap F , k' is the 
momentum measured from the Fermi point. It can be 
written as 



9a W " ^^^ e 



— zfco + avpk' 



(x) 



h (21) 

with r^(x) with the same decay properties as gi^x) 
with an extra factor 2 h ; therefore, the more we are close 
to the Fermi momenta (i.e. —h is large), the more r^(x) 
is a small correction and the propagator is essentially co- 
inciding with the one of a massless Dirac particle. Fi- 
nally V( h > is the effective potential expressed by a sum 
of monomials of order n in the fields i/A-' 1 ) multiplied 
times a kernel wig, while flW is sum of monomials of 

order n in ^ and m in AW with kernels Wi^. Accord- 
ing to power counting, using that <? a (k)^ ~ 2 _,! and 
/ dkgi^k) - 2' 1 , the " naive" scaling dimension of such 
monomials is 



D = 2-n/2 



m 



(22) 



In the RG analysis we have to decompose V^ h ' (a similar 
decomposition must be done also for flw) as 



V {h) = CV {h) + KV {h) 



(23) 



with TZ = 1 — C; CV^ is the relevant or marginal part 
of the effective interaction while TZV^ h ' is the irrelevant 
part. Generally this decomposition is dictated by the 
naive scaling dimension (|22|) ; C should select the terms 
with positive or vanishing dimension D. However, if the 
temperature verifies the condition 



2^ > t h , 



(24) 



where th is the hopping at scale h, there is an improve- 
ment with respect to naive power counting, and certain 
terms which are dimensionally relevant or marginal are 
indeed irrelevant. In order to verify this fact, we can 

where 



split the kernels as W^m = Wn, m 



(b)(h) 



(a)(h) 



is obtained from Wn%, setting t — 0. In the 



(19) case n = 4, m = (with vanishing scaling dimension) 



(a)(h) 



(25) 



so that 



KW™(k') = [Wif h) (k') - Hto )W (0)] + wiTH*) 

(26) 

The first term in the r.h.s. of (1261) can be rewritten as 
k' -dwi^m h \ and this produces an improvement ~ 2 h ~ h 
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in the bound of the kernel, if b! is the scale of the momen- 
tum, which is sufficient to make it irrelevant. Similarly 
the second term in ([26| . namely w\ g (k'), has an ex- 
tra t h 2~ h < 2 h P~ h with respect to the bound for W$, 
which again is enough to make it irrelevant; therefore, 
the true marginal contribution is given by the r.h.s. of 
(|25[) . Therefore the only marginal quartic terms involve 
fcrmions with the same chain index, and that the corre- 
sponding effective coupling coincide with the one of the 
uncoupled t = case. 

Similarly we define, for the terms with n — 2 and the 
same chain index 



-k'8W$ (h) (0) (27) 



Note that W^ h ^ has an extra (th2~ h ) 2 (there are no 
terms linear in th by conservation of the chain index). 
Finally, if n = 2, to = and the fermionic fields have 
different chain index 



(28) 



Note that the terms with n = 2 and an extra deriva- 
tive are irrelevant as they have at least an extra th2T h . 
Therefore 

i— 1 

a— ± a— ± 

a=± J 

In the above expression A^ represents the effective in- 
teraction at momentum scale h, 5h the effective Fermi 
velocity, Uh the shift of the chemical potential and th the 
effective hopping. By definition, Zh,Xh,Vh, &h are exactly 
the same as in the theory with t = 0. It is possible to 
choose v so that i>h remain small for any h. By combining 
Ward-Identities at each Rcnormalization group iteration 
together with Schwinger-Dyson equation it follows, see 
[H, that 



— OO 3— OO (A) 



(30) 



with A_oo(A), (S-oo(A) analytic functions of A; moreover 



Z h ~2 



—r)h 



(31) 



with r\ analytic in A and r\ = aX 2 + 0(A 3 ) with a > 
0. Moreover, in [26| (and references therein) it is 
also proven that kernels Wn h m are analytic functions 
of {X k ,Vk,8 k ,t k } k > h : analyticity (implying the "non- 
perturbative" nature of the method) is a very non triv- 
ial property obtained exploiting anticommutativity prop- 
erties of Grassmann variables, via Gram inequality for 



determinants and Bridges-Battle-Federbush formula for 
truncated expectations. 

Regarding the flow of th we obtain 



th-i 



Z h 



Zh-i 



{t h + A h) ) 



(32) 



with \ j3[ h) \ < Cit h X 2 (t h 2- h ) 2 . It is easy to see by induc- 
tion that \ZhXh — t\ < C2t|A|. Assume indeed that it is 
true for k > h; therefore for A, t small enough 



o 

- i| <2tC 1 X 2 ^ / (t k 2 

k=h 



-k\2 



(33) 



from which, using (|24[) . the inductive assumption follows. 
Note that the effective hopping, even if relevant in the RG 
sense according to naive power counting, remains small 
in this region of temperatures. Moreover, from (1241) we 
obtain the condition between the temperature and the 
hopping 

0~ l >t^(l + 0(X 2 )) (34) 
Regarding the effective source , we define 

CW${ (k',p) = W$$ (h) (0, 0) = 1 (35) 

Indeed the graphs contributing to W^i (0) are one par- 
ticle reducible (as the interaction involves only fermions 
from the same chain) and g^Ck')\k'=o = 0. Therefore 
(assuming that has small support around p = 0) 



(36) 



r 

1212 J M x«i,>x,2, a +<2Al !a ) 



= 1 a=± 



As the flow of the effective parameters corresponding to 
the relevant and marginal operators is bounded, the fol- 
lowing bound is obtained, for /3 verifying (|34[) 



1 



d^Wi% (x)| < C2 h{ - 2 - 



(37) 



In order to compute the conductivity we have to separate 
the terms proportional to t 2 in both the paramagnetic 
and diamagnetic contributions to (|12[) from the rest. We 
write the current-current correlation as 



t- 2 H t (x) 



d 2 W Q 



dA^dA 



X -\o+t- 2 H t (x.) (38) 



where the first term in the r.h.s. is independent from t 
and, from (l37l) 



2»n2 



J dx|£ t (x)| < Ct 2 (|r) 2 ^ 4 ^ 2i 2 C(^ 1 " 

h—ha 



(39) 
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In order to compute the conductivity we still have to com- 
pute < >; introducing the generating functional 

e Wt,x{J) = J p(d^j)e- v ^^ dxJ ^ (40) 
where h x = V>x.i^ x ,2 + ^V'x,! we S et 



< Jx >= n t lo = t I dxi 



9Jx 



<9J x <9J Xl 



|o + A (41) 



where 



n=3 17 17 



dX " 9 Jx^x,..aJxJ° (42) 



and only n odd contribute. From the analogue of 
(|37|) the l.h.s. is bounded by the sum over h of 

^oo =3fi+ l 2 -/ l (, i -l) z -4 g0 that) for t/? smaJ1 







|A| < Cit 2 {^ h tZ- 2 ) 2 < C^t 2 (t^- 2ri ) 2 (43) 



h— hp 

Note finally that 

<9 2 Wo, 



o+ / QXi - ; - T |o (44) 



<9 J x <9 J Xl 



can be rewritten as 

(V>x,l^x,2 - i>i,2<l>Z,2\ V>x,lV>x.2 ~ V , x,2C,2> A + 

or equivalently 

(^x.lV'x.a - V^,2^x,2; ^x,l^x,2 - ^x :2 V'x,2) 0iA ~ 

rf y (V'ilCa - ^x : 2^x,2; ^y,l^y,2 - V'x.aV'y^o.A 

This means that there is an important cancellation be- 
tween the paramagnetic and diamagnetic part of the con- 
ductivity; indeed 

t 2 J dx(e^° - 1) <jx>iy>> M + 0{t 2 {t^) 2 ) 
where 

0x,D;iy,D> 0>A = (C,lV>y,l} 0) A <^y,2^x,2>o,A 

Therefore for oj„ small 

| / dx(e— o - 1) (jx-.DJy^o,, I 



< C\ I dx- ., ., 

'|x|<|«„|-l 1+|X| 2 + 2 '" 



(46) 



rfx 



Ixl^l^l-x 1 + " „ 1 " 



In conclusion the transverse conductivity is given by 
([T]), for t « (i^ 1 « oj n « 1. In the non interacting 
case t~ 2 bj n a^ (w„) ~ — 2 — so that we can conclude that 
the presence of the inter-chain interaction decreases the 
transverse conductivity in this regime. 

The transverse conductivity should be compared with 
the parallel conductivity <rL defined as in (fT2|) with 
J"-",.]' "' replacing' - 
1 one gets 



replacing j D ^,j p ^. For t << f3 1 << u n << 



vfK 



(47) 



where vf = sinpi? + 0(A) is the interacting Fermi ve- 
locity and K is the Luttinger liquid parameter 2rj — 
K + K^ 1 — 2. The above formula, which is part of 
the Haldane's Luttinger liquid conjecture [36j . can be de- 
rived by a Renormalization Group analysis analogue to 
the one described above, see [33] ■ The parallel current- 
current correlation is obtained by a generating func- 
tional similar to (fT3|) in which B(A ± ,iJ;) is replaced by 
/ efoc4 x jx ; after the integration of the fields 
, ip^a^ we get an expression similar to (|20|) 
with B ( - h \A A -, [/Z h ~ip) replaced by B^(A^, y/2%ip) with 



Vi i Via 1 



CB^ h \A\^/z' h ^) = Z i 



(i) 



dxj 



P,U<h) 



and 



Z h 



1 + 0(A) 



(48) 



(49) 



as a consequence of a Ward identity. Therefore the renor- 
malization of the parallel current is proportional to the 
wave function renormalization (while there is no renor- 
malization of the transverse current, see (|35|)) and this 
explains why anomalous power law exponents do not ap- 
pear in the frequency dependence of the parallel conduc- 
tivity. In conclusion, from ([I} and ([4"Tf we see that r\ 
can be independently determined in experiments on two 
chain systems either from the amplitude of the parallel 
conductivity or from the exponent in the orthogonal con- 
ductivity. 

As a final remark, we stress that the above analysis is 

i 

true only for temperatures greater than ~ t 1 -^; at low- 
est temperature the RG analysis would be identical to 
the previous one up to scale t 1 ^ , but at lowest scales 
one should perform the change of variables ([TTf ; the sys- 
tem would be described in terms of two fermions with 
different Fermi momenta (the difference is 0(t T:r *i). In 
this second regime the power counting improvement de- 
scribed above is not valid, and this produces several (not 
a single one, as in spinless Luttinger liquids) effective 
quartic couplings with a generically unbounded RG flow. 
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BILAYER GRAPHENE 

An analysis similar to the previous one can be repeated 
for a model of bilayer graphene, described in terms of 
electrons on the honeycomb lattice interacting through 
an U(l) quantized gauge field, which can represent cither 
the e.m. interaction or the effects of ripples or disorder, 
see e.g. f2l\ . 

We introduce creation and annihilation fermionic op- 
erators tf± = (a% t M +s J = \B\^ dk ^ . ^ 
for electrons with plane index % = 1, 2 and sitting at the 
sites of the two triangular sublattices A^ and As of a 
honeycomb lattice; we assume that A^ has basis vectors 
ha = |(3,±V3) and that A B = A*A + $j, with 5i = (1,0) 
and S2.3 = |(— l,±Vo) the nearest neighbor vectors; B 
is the first Brillouin zone and \B\ = |^=. In the absence 
of e.m. interaction, the Hamiltonian is 

H = H 1 +H 2 + P (50) 

where 

H * = ~ E a ^- +Sj l + c.c. (51) 
■SeA A 

.7=1,2,3 

describes the hopping of fermions in the plane and 

F = - J E [4,i%,2 + 4,2%,i + 

x<EA A 

b + - b~ + b + - b~ - 1 (52) 

describes the fermionic hopping from one plane to an- 
other; either e and t will be assumed small. The inter- 
action with a transverse classical e.m. field is introduced 
via the Peierls substitution. If A is a classical external 
field 



iAi 



b+ - e "Wi&- . +&+ - e" 



? ](53) 



the paramagnetic and diamagnetic part of the transverse 
current are 

,p, x dH(A) , t , _ 

k+5i,1 z+<5i,2 x .2 e,1 x+S u 2 x+Si,l 

. D ,x g 2 ff(A) + _ 

J 2 = „ ,11 = *[ a ?,l%,2 + 

6+ - b~ - +at aZ. +b + - b~ - 

x+5i,l x+8i,2 x ' 2 x > 1 x+<5i,2 



(54) 



and the transverse conductivity is defined as in (|12l) di- 
vided by , the area of the hexagonal cell of the hon- 
eycomb lattice. 



We assume now that the electrons interact through 
an U(l) gauge field; the current-current correlation is 
obtained from the following generating functional 



P(dtf>)P(dA)e 



-V(tp,A)-B(A ± ,ip) 



(55) 



where ip = (a, b) a couple of Grassmann variables 
(with slight abuse of notation, the Grassmann and the 
fermionic operators are denoted with the same sym- 
bol), P(dip) is the fermionic gaussian integration for ip^. 

(i = 1,2 denotes the plane), k = ko,k, fco = + f )> 
with propagator SijgQs.) with 



iko vo£l*(k) 
voQ(k) iko 



(56) 



and v Q = § and Q(k) = § £? =1 e^~ Sl \ The complex 

dispersion relation Q,(k) vanishes only at the two Fermi 
points p F = (^,±^5g) and close to them assumes the 

form of a relativistic dispersion relation fl(pp + k') ~ 
ik[ ± k' 2 . Moreover 

V= -J dxla+ ia b x+Sj ^ e fo g r^+s6 3 )ds + cc] + 
J dxA^a+^a^ + ,<£+«i (57) 



where / dx = J2 x eA A I dx o and Am = { A i, A i 
boson field with propagator 5ijw(p) with 



(Oh 



tu(p) 



yj l^n + C 2 p 2 



(58) 



where x is a cut-off function forbidding momenta either 
too large and smaller than the temperature. Finally the 
source term is given by 



B{A ± ,ip) = J dxA^j 



_L -P,_L 



(59) 



We proceed exactly as in the previous case, writing the 
photon propagator as sum of propagator more and more 
singular in the infrared region, and the fermionic prop- 
agator as a sum of propagators supported close to the 
two Fermi points pp = (^,±^^), labelled by a quasi 
particle index a = ± (labelling the Fermi points) and by 
an integer h < 0: 



«;(p)= X>< A >(p) ( 6 °) 

h=hp 

9 (k)= 9 «(k) + J2 E^(k-p£) 



h—hfs a—± 



with w w (p) supported in 2' 1 " 1 < |p| < 2 h+1 , sup- 
ported on 2' 1 - 1 < |k - pf \ < 2 h+1 and 3 (1) (k) has sup- 
port far from the Fermi points. 
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Assume that we have integrated out the fields 
ijA 1 \..,A^ h+1 \ij/ h+1 \ h>hf) so that 



e W(A ± ) _ p F h (A ± ) 



(61) 



where P(dA^- h ^ ) is the gauge field integration with prop- 



agator 5i_jSfj 



,w 



<i, 



(p), with w^(p) = EL- 



W 



(fc) 



(p), 



while P(dip(- h ) is the integration of the fermionic field 
-0i Q with propagator Si i'6 a a ig§ h (li — Pp) with, if k' = 
k-p^ ) ,g£^ ) (k') = 

_ XfeOO / -ifc w ft (ifci - ak' 2 ) \ h , p {h)\ 

Z h \ «fc(-i*i - afc 2 ) -«fc y ; • 

(62) 

In Eq. (|62l) x^(k') is a cut-off function with support in 

|k'| < 2 h and |i?L' l) (k')| < C|k'| e for some 6 > 0, while 
Z/j and are, respectively, the effective wave function 
renormalization and Fermi velocity on scale h. 

The effective potential + B^ 1 ' expressed by a sum 
of monomials of order n in the fields ip(- h \ m in A^i-ty 
and I in A 1 - , multiplied by kernels ( . According 

to power counting the naive scaling dimension of such 
monomials is 



D = 3-n- 



1 



(63) 



Again there is a dimensional improvement with respect 
to power counting if we are in a range of temperatures 
larger than the hopping, that is 2 h f > th^ where th is the 

hopping at scale h. We can split the kernels as Wjf 1 ^ l — 

<S 3 + <S where is obtained from W$, 

setting t = 0. We define the £ operator in the following 
way 

Note indeed that the extra t h 2~ h < 2 h ^ h in (0) 
is sufficient to make it irrelevant. Regarding the terms 
quadratic in the gauge fields, CW^oip) — Wo°2$\ty + 
pdWQ a)( £\o), where we have used that W^q\o) has 
an extra (2~ h th) 2 with respect to the naive dimension; 
moreover either Wq2o^(0) and dllg 2*0^(0) are vanish- 
ing as consequence of the gauge symmetry, see . Fi- 
nally the terms quadratic in the fermionic variables, if 
they have the same plane index then £W^ (k') = 
^2,o,o°( ) + k'dW^ h \k') where we have used that in 

^2,0,0^ there is an extra gain 0((th2~ h ) 2 , due to the 
conservation of the plane index i. On the other hand for 
the quadratic terms with different plane index 



cww w = w$ (o) 



(ft) 



(65) 



Therefore 



CVW(A,if>)=t h I dxjf- ± + 
dk dp 



(66) 



B\ (27T) 



where: e ,h = e ,ft, &i,h = v h e ith , = e 2 ,h (thanks 
to discrete rotational symmetry), := T fl (^,0) (with 
rjj" = -il, TY = -U2, = -cjo-i and o- 1)2 the first 
two Pauli matrices) and IZV^ a sum of terms that are 
irrelevant in the RG sense. By construction the flow of 
the effective parameters is the same as in the model with 
t = 0; it was show in [3(|, by a rigorous implementation 
of Ward Identities in the RG scheme, that the effective 
charges flows to a line of fixed points and the Fermi ve- 
locity increases up to the light velocity 



e h -> e_ 



(67) 



Moreover, the wave function renormalization Zh diverges 
with anomalous exponents 



12tt 2 



Finally regarding the flow of th we obtain 



1-/1-1 



Zh-l 



{t h + A h) ) 



(68) 



(69) 



with |/3 t ^ | < Cie 6 t/j[|fr] 2 , and again by induction \ZhXh~ 
t\ < C^te 6 . We assume that the temperature verifies 
(gH which implies /3" 1 > t^n (1 + 0(e 2 )). Regard- 
ing the effective source B^ h \ we define CW^q 1 (k', p) = 
W201 (0,0) = 1 as again the graphs contributing to 

W 2 *o 1*^(0) are one particle reducible and (k')|k'=o = 
0. As the flow of the effective parameters corresponding 
to the relevant and marginal operators is bounded, the 
following bound is obtained, for h > hp (order by order 
in the renormalized expansion) 



(70) 



Using the same notation as in ([38 



J dx|so||# t (x)| < Ct 2 (|r) 2 ^ 4 < 2t 2 C(t/3 1 - 2 ") 2 

h=hp 

(71) 

Moreover, as in the previous case we introduce a gener- 
ating functional VVt. e (J) with source t J dxJ x h x where 
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From the analogue of (|70[) the l.h.s. is bounded by the 
sum over h of £^ =3 p+^-M*- 2 )^ 4 so that, for t/3 
small 



X |A| < t 2 J2 /32 h [^V] 2 < Ct 2 (tj3 1 ~ 2r) ) 2 (73) 



-/i 



h—ha 



z 



Note finally that 

\Jx >Jy /o, e V* 'Jy /o,e 

and 



(74) 



dyi\x \ 



FO^n 



|x|<fj~ 



1 + |x| 4 + 2 " 

1 



Cl W d3xKI UN^- „ 

Therefore the conductivity in the interacting case is given 
by © for t « « cj n « 1, that is the trans- 
verse conductivity decreases with the frequency with the 
anomalous exponent 2rj. In absence of planar interac- 
tion t~ 2 (jp(u n ) ~ i, so that we can conclude that the 
presence of planar long range interaction producing Lut- 
tinger liquid behavior decreases the transverse conduc- 
tivity. Note also that the parallel conductivity does not 
display any anomalous power law, as a consequence of a 
Ward Identity implying the analogue of (|49|) . see [30j. 



APPENDIX: THE NON INTERACTING CASE 



In the case of the two chain model if A = 



t 2 LJ n <Tp(u} n ) 



dkg(k + p) 5 (k) 



dkg(k)g(k)}\ p=0 + 0{{pt) 2 ) 



Note that 



«l im 7^I>(k + p)ff(k)| p= o =0 

p,L— >oo Lip e —^~ 



while 



(75) 



(76) 



In the case of bilayer graphene, we get t 2 a^(iu n ) — 
— [ ^dk[F(k,k + p)-F(k,k)}\p=o + 0((tf3) 2 ) 

UJ n J (^7Tj 



where 



F(ki,k 2 ) = 2[pxi(ki)5ii(k 2 ) + 522(ki)3 2 2(k 2 ) 
+5i 2 (ki)5i 2 (k 2 ) + .92i(ki)g2i(k 2 )] (77) 



The first term in the r.h.s. can be written as its value 
for (3 = oo plus a rest 0(/3 _1 ); the integral in the limit 
/3 = oo can be decomposed in a part integrated in the 
region < e and |^(fc)| > e; the second term is 

vanishing for u) n — while in the first the contribution 
from the first two terms in (1771) are vanishing by parity, 
while the rest gives ^ at vanishing external frequency. 
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